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We prove a theorem about magnification relations for all generic general caustic singularities up to 
codimension five: folds, cusps, swallowtail, elliptic umbilic, hyperbolic umbilic, butterfly, parabolic 
umbilic, wigwam, symbolic umbilic, 2°^^ elliptic umbilic, and 2""^ hyperbolic umbilic. Specifically, 
, . . we prove that for a generic family of general mappings between planes exhibiting any of these 

0^ I singularities, and for a point in the target lying anywhere in the region giving rise to the maximum 

. number of real pre-images (lensed images), the total signed magnification of the pre- images will 

' always sum to zero. The proof is algebraic in nature and makes repeated use of the Euler trace 

04 ' formula. We also prove a general algebraic result about polynomials, which we show yields an 

{^jQ, interesting corollary about Newton sums that in turn readily implies the Euler trace formula. The 

■ wide field imaging surveys slated to be conducted by the Large Synoptic Survey Telescope are 

expected to find observational evidence for many of these higher-order caustic singularities. Finally, 
since the results of the paper are for generic general mappings, not just generic lensing maps, the 
^SJ ■ findings are expected to be applicable not only to gravitational lensing, but to any system in which 

these singularities appear. 
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I. INTRODUCTION 



, One of the key signatures of gravitational lensing is the occurrence of multiple images of lensed sources. The 
^ • magnifications of the images in turn are also known to obey certain relations. These relations fall into two types: 
\^ [ "global" and "local". "Global" magnification relations involve all the images of a given source, but they are not 
CO . universal because the relations depend on the specific class of lens models used. Examples of such relations can be 
CN ■ found in Fetters et al. 2001 [Ml, p. 191], Witt & Mao 1995 [H, Rhie 1997 (H], Dalai 1998 'sl, Witt & Mao 2000 

! d^, Dalai & Rabin 2001 and Hunter & Evans 2001 [l2|. As shown in Werner 2007 such relations are in fact 
Tj" ■ topological invariants. 

' By contrast, "local" magnification relations are universal, but they apply only to a subset of the total number 
' of images produced. Two well-known examples of local magnification relations are the fold and cusp relations. For 
a source near a fold or cusp caustic, the resulting images close to the critical curve are close doublets and triplets 
whose signed magnifications always sum to zero (e.g., Blandford & Narayan 1986 [6|, Schneider & Weiss 1992 [25j . 
Zakharov 1999 [30|, [Ml, Chap. 9]). The universality of these relations means that they hold independently of the 
choice of lens model. In addition, the fold and cusp relations have been shown to provide powerful diagnostic tools 
' for detecting dark substructure on galactic scales using quadruple lensed images of quasars (e.g., Mao & Schneider 
■ 1998 [3, Keeton, Gaudi & Fetters 2003 and 2005 [ll,ll|). 

Recently, Aazami & Fetters 2009 [1], which we consider Faper I, established a universal magnification theorem 
for some of the higher-order caustics beyond folds and cusps, namely, the swallowtail, elliptic umbilic, and hyperbolic 
umbilic singularities. These are generic caustic surfaces or big caustics occurring in a three-parameter space. Slices 
of the big caustics give rise to generic caustic metamorphoses (e.g., [M|, Chapters 7 and 9), all of which occur in 
gravitational lensing (e.g., Blandford 1990 "5], Fetters 1993 [20|, Schneider, Ehlers, & Falco 1992 and [MI)- 
It was shown in [l[ that for lensing maps close to elliptic umbilic and hyperbolic umbilic caustics, and for general 
mappings exhibiting swallowtail, elliptic umbilic, and hyperbolic umbilic caustics, the total signed magnification for 
a source lying anywhere in the region giving rise to the maximum number of lensed images, is identically zero. As an 
application, they used the hyperbolic umbilic to show how such magnification relations may be used for substructure 
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studies of four-image lens galaxies. 

The proof of these relations in [l| was elementary, but long, and thus was not amenable to higher-order caustics 
beyond the three mentioned above. An elegant geometric technique, based on Lefschetz fixed point theory, has since 



been employed on these three singularities by Werner 2009 [27[. The aim of our current paper is to extend these 
results to all the remaining higher-order generic general singularities up to codimension 5. In other words, our findings 
are expected to be applicable not only to gravitational lensing, but to any system where these singularities appear. 

We prove that to each generic general caustic singularity of codimension up to 5 — not just the fold, cusp, swallowtail, 
elliptic umbilic, and hyperbolic umbilic, but also the butterfly, parabolic umbilic, wigwam, symbolic umbilic, 2"'' elliptic 
umbilic, and 2"'' hyperbolic umbilic — is associated a magnification sum relation of the form 

^ on. = . 

In other words, for generic families of general mappings between planes exhibiting these singularities, and for a point 
anywhere in the region of the target space giving rise to the maximum number of lensed images, the total signed 
magnification is identically zero. Furthermore, as emphasized in magnification sum relations are in fact geometric 
invariants, because they are the reciprocals of Gaussian curvatures at critical points. 

Shin & Evans 2007 p/Fj] constructed a realistic lens model for the Milky Way Galaxy and showed that it exhibited 
butterfiy caustics (see also Evans & Witt llSll for another class of lens models that exhibit butterfiy caustics) . More 
recently, Orban de Xivry & Marshall 2009 [1^ created an atlas of predicted gravitational lensing due to galaxies and 
clusters of galaxies that can exhibit several of these higher-order caustic singularities, and estimated the probabilities 
for their occurrence. They showed how a galaxy lens with a misaligned disk and bulge can generate swallowtails and 
butterfiies, two merging galaxies or galactic binaries can produce elliptic umbilics, and galactic clusters can create 
hyperbolic umbilics. These lensing effects are expected to be seen by the Large Synoptic Survey Telescope. 

Concerning the tools of the paper, we mention that the Euler trace formula was employed in ^ to determine "global" 
magnification relations for special classes of lens models. They used an analytical approach whereby they derived the 
Euler trace formula using residue calculus. We show that the Euler trace formula also lends itself quite naturally to 
"local" magnification relations for generic general caustics, not just those occurring in gravitational lensing. In fact, 
along with our main theorem, we also prove a general algebraic theorem about polynomials (Proposition [2]) and derive 
a result about Newton sums as a corollary, which we show implies the Euler trace formula. Our main theorem is not 
a direct consequence of the Euler-Jacobi formula, of multi-dimensional residue integral methods, or of Lefschetz fixed 
point theory, because some of the singularities have fixed points at infinity. 

The outline of the paper is as follows. Section [III reviews the necessary singular-theoretic terminologies and results. 
Section Hill states our main theorem, which is for general mappings. As preparation for the proof of our main theorem, 
in Section IIVI we establish a recurrence relation for the coefficients of the unique polynomials in cosets of certain 
quotient rings and show that this relation yields a fact about Newton sums which can be employed to readily obtain 
the Euler trace formula. We then use the results of Section |IV] to prove the main theorem in Section |Vl 



II. HIGHER-ORDER CAUSTIC SINGULARITIES 



In what follows, the term "universal" or "generic" is used to denote a property that holds for an open, dense 
subset of mappings in the given space of mappings. With that said, consider a smooth, n-parameter family ^^^(x) 
of functions on an open subset of that induces a smooth (n — 2)-parameter family of mappings fc(x) between 
planes (n > 2). Given fc(x) — s, call x the pre-image of the target point s. Critical points of fc are those x for 
which det(Jacfc)(x) = 0. Generically, the locus of critical points forms curves called critical curves. The value fc(x) 
of a critical point x under fc is called a caustic point. These typically form curves, but could be isolated points. 
Varying c causes the caustic curves to evolve with c. This traces out a caustic surface, called a big caustic, in the 
n-dimensional space {c, s} — R"^^ x M^. Beyond the familiar folds and cusps, these surfaces form higher-order caustics 
that are classified into universal or generic types for locally stable families fc. Generic c-slices of these big caustics 
are commonly called caustic metamorphoses. 

The universal form of the {n — 2)-parameter family fc is obtained by using F^-^s to construct catastrophe manifolds 
that are projected into the space {c, s} = K"~^ x to obtain local coordinates for fc (e.g., Majthay 1985 \T^ . 
Castrigiano & Hayes 1993 7], Golubitsky & Guillemin 1973 [HI). These projections of the catastrophe manifolds are 
called catastrophe maps or Lagrangian maps, and they are differentiably equivalent to fc (see (2ll . pp. 273-275]). The 
locally stable families Fc,s and their induced maps fc are generic for n < 5, and have caustic singularities that are 
classified according to the parameter n. For n = 3, the singularities generically divide into three types: swallowtails, 
elliptic umbilics, and hyperbolic umbilics. When n — A, the singularities generically divide into two types: butterflies 
and parabolic umbilics. For n = 5, they divide into four types: wigwams, symbolic umbilics, 2"'' elliptic umbilics, and 
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2"'' hyperbolic umbilics; see Table ID A detailed treatment of these issues can be found in Poston and Stewart 1978 
[H, Gilmore 1981 [H, [H, Arnold 1986 Q, and HH, Chap. 7]. 

III. STATEMENT OF MAIN THEOREM 

For the generic families of functions i^c,s in Table HJ we define the magnification 97l(xi;s) at a critical point x,; of 
fc,s relative to x = (x, y) by the reciprocal of the Gaussian curvature at the point (x^, _Fc,s(xi)) in the graph of Fc-y. 

^{^^■:S) = . . • 

Gauss(xi,Fc,s(xi)) 

This makes it clear that the magnification invariants established in our theorem are geometric invariants. 

Theorem 1. For any of the universal, smooth n-parameter family of general functions F^s {or general mappings f^) 
in Table\^ and for any non-caustic point s (light source position) in the indicated region, the following results hold 
for m^ = m{x.r,s): 

1. A2 (Fold) Magnification relations in two-image region: 

97li + 97I2 = . 

2. A3 (Cusp) Magnification relations in three-image region: 

Til TI2 + TI3 ^ . 

3. A4 (Swallowtail) Magnification relation in four-image region: 

SWi + 97I2 + SW3 + 9714 = . 

4- (Elliptic Umbilic) Magnification relations in four-image region: 

OTi + 97l2 + OTa + 97I4 = . 

5. (Hyperbolic Umbilic) Magnification relations in four-image region: 

Srjli + 97I2 + 97l3 + 2^4 = . 

6. A5 (Butterfly) Magnification relation in five-image region: 

Tli+m2+ TI3 + 9714 + 9JT5 = . 

7. (Parabolic Umbilic) Magnification relations in five-image region: 

OTi + 9Jl2 + 9JI3 + 9JI4 + 9JI5 = . 

8. Ag (Wigwam) Magnification relations in six-image region: 

OTi + 97t2 + 9Jl3 + M4 + aJls + OTe = . 

9. Eq (Symbolic Umbilic) Magnification relation in six-image region: 

9Jli + 97l2 + 97I3 + 97l4 + Sr«5 + aJlg = . 

10. Dq (2"'^ Elliptic Umbilic) Magnification relations in six-image region: 

mil + 9Jt2 + 9Jl3 + 9714 + 97l5 + 2)^6 = . 
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Fold (2D) 


Fs{x,y) = -six + S2y+ \x'^ - \y^ 
i{x,y) = {x , 2/2) 


Cusp (2D) 


Fs{x, y) = -Six + S2y + \x'^ - 5S12/2 - 
i{x,y) = {x , xy + y'^) 


Elliptic Umbilic (3D) 


Fc,six,y) = Six + S2y + cix^ + 2/2) + x^ - "ixy^ 
^c(x, y) = (3j/2 - 3a;2 _ 2cx , 6xy - 2cy) 


Hyperbolic Umbilic (3D) 


Fc,a{x, y) = Six + S2y + cxy + x'^ + y'^ 
^c(x,y) = {-3x^ - cy , -Zy^ - cx) 


Swallowtail (3D) 


Fc,s(x,y) = Six - S2y - ^S2X^ + ^y^ - ^cx^ - ^x^ 
^c{x,y) = (xy + ca;2 + X* , y) 


Butterfly (4D) 


Fc,s{x, y) = x^ + cix* + C2X^ + S2x'^ + six + \y'^ — S2y 
fc(x, y) = (—2x2/ — 3c2x2 — 4cix^ — 6x^ , y) 


Parabolic umbilic (4D) 


Fc,s(x,y) = x^y + y* + cix^ + 022/^ - six - S2y 
fc(x, y) = (2cix + 2xy , 2c2y + + 4y^) 


Wigwam (5D) 


Fc,six, y) = x"^ + ClX'' + C2X'^ + C3X-^ + S2X'^ + Six + _ g^y 

fc(x,2/) = (-2x2/ - 3c3x2 - 4c2X^ - Scix** - 7x® , y) 


Symbolic umbilic (5D) 


Fc,s{x, 2/) = x^ + 2/'* + C1X2/2 + C2xy + C3y'^ + S22/ + six 
fc(x,2/) = (-3x2 _ ^^y2 _ ^ _^y3 _ 2ci Xy - C2X - 2c32/) 


2"^^ Elliptic umbilic (5D) 


Fc,s{x, y) = x^y -y^ + ciy^ + C2y^ + csy"^ + S2y + six 
fc(x, y) = (-2x2/ ) -x^ + 5y* - 4ciy^ - 3c22/^ — 2c32/) 


2"d Hyperbolic umbilic (5D) 


Fc,s(.x, y) = x'^y + y^ + cij/* + C22/^ + 032/^ + S22/ + six 
fc(x, y) = (-2x2/ , -x'^ - 52/" - 4cij/3 - 3c22/2 - 2c32/) 



TABLE I: For each type of caustic singularity listed, the second column shows the universal local forms of the smooth n- 
parameter family of general functions -Fes, along with their (n — 2)-parameter family of induced general maps fc between 
planes. The numbers 2D, ZD, etc., denote the codimension of the given singularity. 
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11. Dq (2"^^ Hyperbolic Umbilic) Magnification relations in six-image region: 

Ml + OTa + + 97^4 + 2^5 + a^^e = . 

We use the A, D, E classification notation of Arnold 1973 in the theorem. This notation highlights a deep link 
between the above singularities and Coxeter-Dynkin diagrams appearing in the theory of simple Lie algebras. As 
mentioned in the introduction, the fold and cusp magnification relations are known [i,[2l|,[2i,[3q]. The magnification 
relations for the swallowtail, elliptic umbilic, and hyperbolic umbilic were discovered recently in 

Remark. The results of Theorem [T] actually apply even when the non-caustic point s is not in the maximum number 
of pre-images region. However, complex pre-images will appear, which are unphysical in gravitational lensing. 



IV. A RECURSIVE RELATION FOR COEFFICIENTS OF COSET POLYNOMIALS 

In this section, we present some notation and a proposition about polynomials that will yield the Euler trace formula 
as a corollary. The notation and the latter are used in the proof of Theorem [TJ The proposition itself is proved in the 
Appendix. 

We begin with some notation. Let C[a;] be the ring of polynomials over C and consider a polynomial 

^p{x) = a„a;" + • ■ • + aix + ag G C[x] . 

Suppose that the n zeros Xi, . . . ,Xn pf (p{x) are distinct (generically, the roots of a polynomial are distinct) and let 
ip'{x) be the derivative of ip{x). Also, let R C C(a;) denote the subring of rational functions that are defined at the 
roots Xi of f{x): 

I pixj 

R — s , : p{x), q{x) e C[a;] and q{xi) ^ for all roots Xi 

Let (1^3(2;)) be the ideal in R generated by ip{x) and denote the cosets of the quotient ring R/(ip{x)) using an overbar. 
Below are two basic results that we prove in the Appendix (see Claim 2) for the convenience of the reader: 

• Members of the same coset in R/{(p{x)) agree on the roots Xi of (p{x), that is, if hi{x) and h2{x) belong to the 
same coset, then hi{xi) = h2{xi). 



• Every rational function h{x) e R has in its coset h{x) £ R/{ip{x)) a unique polynomial representative of 
degree less than n. 

Proposition 2. Consider any polynomial (p{x) = a„a;" + ■ • ■ + aix + ao G C[x] with distinct roots and any rational 
function h{x) G R. Let 

h^,{x) — Cn-ix"^^ + • • • + CiX + Co 



be the unique polynomial representative of the coset h{x) G R/{ip{x)) and let 

r{x) = bn-ix''^^ H \-bix + bo 



be the unique polynomial representative of the coset (p'{x)h{x) G R/{ip{x)). Then the coefficients of rix) are given in 
terms of the coefficients of h^:{x) and ^p{x) through the following recursive relation: 

bn-i = Cn-ibn-'Ln-l -\ + Cl6„-i,l + Co6n-i,0 i = l,...,n, (1) 

with 

bn-i,o = [n - (* - 1)) an-(i-i) , i = 1, . . . , n , 

bn-i,k = bn-l.k-l + bn-(i+l).k~l , i = 1, . . . ,71 , fc=l,...,n-l, 

an 

where b-i^k-i = 0. 



(2) 
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By Proposition [21 if rfe(x) is the unique polynomial representative of the coset (f'{x)x^ S R/{(p{x)), then 

rfc(a;) = 6„_i,fex""^ H \-biMX + bQ^k, fc = 0, 1, . . . , rt - 1 , (3) 

where its coefficients are given in terms of the coefficients of (p{x) through 

Corollary 3. Assume the hypotheses and notation of Proposition\^ Given the distinct roots Xi, . . . , a;„ of (p{x), the 
Newton sums Nk = satisfy: 



fc = 0, 1, . . . , 71 — 1 



(4) 



In other words, the quantity a„ Nk equals the {n — l)st coefficient of the unique polynomial representative ([3]) of the 
coset (/?'(x)a;*-' in R/{ip{x)). 

Proof. Note that for k — 0, eqn. ([2]) in Proposition [2] yields 

For 1 < fc < n — 1, there is a known recursive relation for Nk, in terms of A^i, N2, ■ ■ ■ , Nk-i; see, e.g., Barbeau 1989 [3, 
p. 203]. It is given by 



kan^k + On-k+lNi + an-k+2N2 H h a„_lA^fc-l + OnNk = 



(5) 



We proceed by induction on fc for 1 < /c < n — 1. For k — 1, eqn. ([5]) implies A^i = — °^ ^ , while eqn. ^ gives 
bn-1,1 = —On-i = OnNi, which agrees with eqn. (jj]). Now assume that bn-ij — OnNj for j — 1, . . . , k — 1. To 
establish the result for j — k, we shall repeatedly apply Proposition [51 



bn-i,k — — 



an 



bn-l,k-l + &n-2,fe-l 

, , ari-2 

On-l,fc-l 



bn-l.k-2 + bn-3,k-2 



a, 



On-l,k-l On-l,k-2 + 



n 



bn-l,k-3 + b„-4 k-3 



flri-l , 0,11-2 , 0,n-3 , lT!-(fc-l) , Un-k , , 

— On-l,k-l 0„-l,fc-2 On-lM-3 — ■ ■ • On-1,1 On^lfi + On-(fc+l),0 • 

Oji Ojyi tin On Ojii 

= - {a-n-lNk-l + an-2Nk-2 + an-3Nk-3 H \- an-{k-l)Ni + kOn-k) 

= anNk , 

where 6n-i,o = non and &„-(fc+i) = [n — k)an-k follow from eqn. ([2]) in Proposition [21 and the last equality is due 
to ®. □ 

Corollary 4 (Euler Trace Formula). Assume the hypotheses and notation of Proposition\^ For any rational function 
h(x) e R, the following holds: 



^h{xi) 



bn-l 



(6) 



where is the (n — l)st coefficient of the unique polynomial representative r(x) of the coset ip'{x) h(x) G R/ {(f{x)) 

and On the nth coefficient of (p(x). 



Proof. Let h^,{x) be the unique polynomial representative of the coset h{x) G R/{(p{x)). First note that, since h{x) 
and h^(x) belong to the same coset, we have h(xi) — h^[xi). The Euler trace formula now proceeds from a simple 
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application of Propositon[2]and Corollary [3j 

n n n n— 1 n — 1 n n— 1 

i=l i=l 2=1 j=0 j=0 i=l j=0 

= c„_i7V„_i + • • • + ciiVi + coiVo 

= c„_i f ^lizi-^) + ...+C, f ^) + CO f ^) (by Corollary© 

\ Cin J \ an J \ an J 

_ Cn — 1 — 1 , n — 1 ~t~ 

an 

bn-1 



(by Proposition [2) □ 



Remark. Dalai & Rabin 2001 Q gave a different proof of the Euler trace formula, one employing residues. 

V. PROOF OF THE MAIN THEOREM 

We begin by establishing some preliminaries before starting the computational part of the proof. Given a family of 
functions F^.s, a parameter vector (co,So) is called a caustic point of the family if there is at least one critical point 
Xq of fco,so (i-e., xq satifies grad-F'co.so(xo) = 0) such that the Gaussian curvature at (xq, -Fc,s(xo)) in the graph of 
Fc^s vanishes. Furthermore, for the list of singularities in Table HI the mappings fc are induced by the families -Fes- 
In fact, a direct computation shows that for all the singularities, we can obtain fc through the gradient of i^c,s as 
follows: 

gradFc,s(x) = fc(x) = s . 

We can also express the magnification in terms of the general mappings fc induced from the n-parameter family of 
functions -Fc.s- To do so, recall that the Gaussian curvature at a point (x, Fc.s(x)) in the graph of Fc.s is given by 

Gauss(x,Fc,(x))= dot (Hess Fc,)(x) 



l + |gradFc,s(x)|2 
At a critical point xq , the magnification of xo is then given by 

1 1 



OT(xo;s) = 



Gauss(xo, Fc,s(xo)) det(Hess i^c,s)(xo) 



Note that caustics are characterized by the family F^- s having at least one infinitely magnified critical point. A 
computation also shows that for all the singularities in Table HI the following holds: 

det(Jacfc) = det(HessFc,s) • 
Consequently, we can also express the magnification at a pre-image Xj of s under fc as 

m = OJl(x,;s) = ■ ^„ - , fc(x,) - s . 

det(Jacfc)(xj) 

Observe that caustics in the target plane of fc are given equivalently as points s where the Jacobian determinant of 
fc vanishes. Now, given an induced mapping fc and a target point s = (si, S2), we can use the pair of equations 

(si,S2) = ic{x,y) = {fi^c{x,y),f2A^^y)) 

to solve for (x, y) in terms of (si, S2), which will give the pre-images x^ — (xi, yi) of s under fc. 

For the singularities in Table [H we shall see that the pre-images can be determined from solutions of a polynomial 
in one variable, which is obtained by eliminating one of the pre-image coordinates, say y. In doing so we obtain a 
polynomial <^(a;) G C[a;] whose roots will be the x-coordinates Xi of the different pre-images under fc: 

ip{x) = Gnx"' + ■ ■ ■ + aix + ao . 
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Generically, we can assume that the roots of (p{x) are distinct, an assumption made throughout the paper. 

We would then be able to express the magnification dJl{x, y; s) at a general pre-image point [x, y) as a function of 
one variable, in this case x, so that 

9Jl(x, y{x)]s) = = = m{x) , 

J{x,y{x)) J{x) 

where J = det(Jacfc) and the explicit notational dependence on s is dropped for simplicity. Since we shall consider 
only non-caustic target points s giving rise to pre-images {xi,y{xi)), we know that J{xi) ^ 0. Furthermore, we shall 
only consider non-caustic points that yield the maximum number of pre-images. In addition, for the singularities 
in Table m the rational function 9Jt(x) is defined at the roots of f{x)^ i.e., 9Jl(x) e R. Now, denote by m(x) the 
unique polynomial representative in the coset ^'{x) 97l(x) G Rj {(p{x)), and let &„_i be its (n — l)st coefficient. In the 
notation of Proposition [21 we have h{x) = dJl{x) and r{x) = m.{x). Euler's trace formula (Corollary [H]) then tells us 
immediately that the total signed magnification satisfies 



Vot. = ^. (7) 
— ^ n„ 



It therefore remains to determine the coefficient for each caustic singularity in Table [H Next to each singularity 
below we indicate the value of n — 1, which is the codimension of the singularity. 

Finally, we mention that the full theorem is not a direct consequence of the Euler-Jacobi formula, of multi- 
dimensional residue integral methods, or of Lefschetz fixed point theory, because some of the singularities have 
fixed points at infinity. 

1. Fold (1): Its corresponding induced map f is 

f{x,y) = {x, y^) . 

Let s = (si, S2) be a non-caustic point and let us determine the maximum number of images of s. Setting 
we obtain x = si and find that the y-coordinates of the pre-images are the two real zeros of the polynomial 

¥'(2/) = r/ - S2 ■ 

Consequently, there is a maximum of two pre-images. The magnification, expressed in the one variable y, is 
given by = 1/J{y) — 1/2?/. Since (p'{y) — 2y — J{y), we have 

ip'iy)miy) = Jiy)miy) = 1 . 



But this implies that the unique polynomial representative in the coset ip'{y)'^{y) is the polynomial tn(y) = 1. 
Since the [n — l)st coefficient is 61 = 0, we conclude via ^ that the total signed magnification in the two-image 
region is zero: 

2 , 

^ ^ /In 



2. Cusp (2): Its corresponding induced map f is 

f{x,y) = [x, xy + y^) . 

As with the fold, let s = (si, S2) be a non-caustic point and let us determine the maximum number of images 
of s. Setting 

f{x,y) = (S1,S2) , 

we obtain x — si and find that the y-coordinates of the pre-images are the three real zeros of the polynomial 



=y^ + siy - S2 
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So, there is a maximum of three prc-images. The magnification, expressed in the one variable y, is given by 
Tl{y) = l/J{y) = l/(3y^ + Si). Once again we have (p'{y) = J{y), so that 



As with the fold, it follows that m(y) = 1 and 6„_i = 62 = 0, so that the total signed magnification in the 
three-image region is zero: 



3. Swallowtail (3): Its corresponding 1-parameter induced map fc is 

fc{x, y) = {xy + cx^ + , y) . 
Let s = (si, S2) be a non-caustic point and let us determine the maximum number of images of s. Setting 

f(a:,?/) = (si,S2) , 

we obtain y = S2 and find that the x-coordinates of the pre-images are the four real zeros of the polynomial 



which gives a maximum of four pre-images. The magnification is £!Jl(a;) = 1/J{x) = l/(4a;^ + 2cx + S2) and 
<f'{x) = J{x), so that 

ip'{x)m{x) = J{x)mi{x) = 1 . 
We thus have m(a;) = 1, = 63 = 0, and the total signed magnification in the four-image region is zero: 

4. Elliptic Umbilic (3): Its corresponding 1-parameter induced map fc is 

fc(x, y) = {3y^ — — 2cx , 6xy ~ 2cy) . 

Setting i{x,y) = (si,S2) for a non-caustic point (si,S2), we get y = S2/{Qx — 2c), which we use to get the 
following degree 4 polynomial for the x-coordinates of the pre-images: 

fix) = Ac^si - 3sl + (8c^ - 24csi)x - (SGc^ - 3Qsi)x^ + 108x^ . 

Hence there is a maximum of four pre-images. The magnification is Tl{x,y) = l/(4c^ — 36x^ — 36?/^), which 
becomes a function of x: 

^^""^ " lix) " 4c2 - 12si - 24cx - 72x2 ' 
This time <^'(x) 7^ J{x), but the situation is remedied if we multiply through by 2c — 6x to get 

ip'{x)m{x) = [J{x)mi{x)] (2c - 6x) = 2c - 6x . 

Thus the unique polynomial representative in the coset ip'(x)V}l(x) is the polynomial tn(x) = —6x + 2c. Since 
bn-i = 63 = 0, we have 



5. Hyperbolic Umbilic (3): Its corresponding 1-parameter induced map fg, for a given target point (si, S2) lying 
in the four-image region, is 



f'{y)dniy) = J{y)m{y) = 1 . 




ip{x) = X^ + ex? -\- S2X — Si 



Tii + m2 + m3 + m4 = o . 



fc{x, y) = (-3x^ - cy , -3y^ - cx) = (si, S2) • 
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We eliminate y to obtain a polynomial in the variable a;, given by 

<fi{x) = -3s? - c^S2 - c^x - 18sia;2 - 27a;^. 
The magnification is 9Jt(a;,t/) = !/{—(? + 36a;y). Substituting for y via fc(a;,y) = (si, S2), we obtain 

1 c 
" " -c3 - msix - 108a;3 ' 

It follows that cJ{x) = (p'{x), so that 

ip'{x)m{x) = [J(x)9Jl(x)] c = c . 



Thus the unique polynomial representative in the coset ip'{x)dJl{x) is the polynomial m{x) = c. Since 6„_i = 
63 = 0, we have 

SDTi + mt2 + 9^3 + 3^4 = . 

6. Butterfly (4): Its corresponding 2-parameter induced map fc, for a given target point (si,S2) lying in the 
five-image region, is 

fc(x, y) = {-2xy - Zcix^ - ^c\x^ - , y) = (si, S2) . 
We eliminate y and are left with a polynomial in the variable x, 

ip{x) = —Si — 2s2X — 3c2X^ — 4cix^ — 6x^ , 
whose roots are the x-coordinates of the five pre- images. In this case J{x) = (p'{x), so that 

ip'{x)m{x) = 1 . 



Thus m{x) = 1 is the unique polynomial representative in the coset ^p'{x)^M{x). Since 6„_i = 64 = 0, it follows 
that 

2Jli + OJI2 + 9K3 + 2JI4 + 2Jl5 = . 

7. Parabolic Umbilic (4): Its corresponding 2-parameter induced map fc, for a given target point (si,S2) lying 
in the five-image region, is 

fc{x, y) = {2cix + 2xy , 2c2y + x"^ + V) = (si, S2) . 

We eliminate x to obtain a polynomial in the variable y, 

ip{y) = -sf + 4c?S2 - (8c?C2 - 8ciS2)y - (I6C1C2 - 4s2)y' - (16c? + 8c2)y^ - 32ciy^ - Wy^ . (8) 

The magnification is Tl{x,y) = l/(4ciC2 — + 4c2y + 24ciy^ + 24j/'^). Substituting for x"^ via the equations 
fc{x,y) = (si,S2), we obtain 

m{y) ^ ^ 



J{y) 4ciC2 - 4s2 + 12c2t/ -|- 24ciy2 -|- 40y3 
Although ip'{y) ^ J{y), the situation is remedied if we multiply through by — 2ci — 2y : 

ip'{y)m(y) = [Jiy)miy)] (-2ci - 2y) = -2ci - 2y , (9) 

from which we immediately conclude that m(y) = — 2ci — 2y is the unique polynomial representative in the coset 
if' {y)9Jl{y) . Since bn-i = 64 = 0, we have 

OTi + 9Jt2 + SIJ13 + SIK4 + = 0. 
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8. Wigwam (5): Its corresponding 3-paranieter induced map fc, for a given target point (si,S2) lying in the 
six-image region, is 

fc(x, y) = {-2xy - 3c3X^ - 4c2X^ - 5cix'^ — 7x^ , y) = (si, S2) • 
We eliminate y and are left with a polynomial in the variable x, 

ip{x) = —2s2X — Scaa;^ — 4c2a;^ — Scia;'* — — si , 
whose roots are the a;-coordinates of the six pre- images. In this case J{x) = ip'{x), so that 

ip'{x)m{x) = 1 . 

Therefore, as with the fold, ctisj). s\\ allowtail, and butterfly, we conclude that m(a;) = 1 is the unique polynomial 
representative in the coset if' {x)f}Jl{x) . Since 6„_i = 65 = 0, it follows that 

9Jli + £Crt2 + aJts + 9JI4 + 9Jl5 + SDTe = . 

9. Symbolic Umbilic (5): Its corresponding 3-parameter induced map fc, for a given target point (si,S2) lying 
in the six-image region, is 

fc(a;, y) = {-3x^ - cit/^ - C2y , -4y^ - 2cixy - C2X - 2czy) = (si, S2) . 

We eliminate x to obtain a polynomial in y, 

V{y) = -clsl - 3sl - cly - ■iciC2Siy - 12c3S2y - Bcicly^ - Ucly"^ - Aclsiy"^ 
- 8c?C2y^ - 24s2y^ - 4c?/ - 48C3/ - 48/ . 

The magnification is DJl{x,y) = l/{—c^ + 12c3X + 12c\x'^ — Ac\C2y — ^c^y^ + 72a;/). Substituting for x via the 
equations fc(a;, y) = (si,S2) gives 



1 C2 + 2ciy 



J{y) -4 - 4C1C2S2 - 12C3S2 - (lOcici + 24c§ + 8c?S2)j/ - (24cfc2 + 72s2)/ - (16c? + 192c3)/ - 288/ 
It is not difiicult to check that ifi'{y) = J{y){2ciy -h C2), so that 

ip'{y)Tl{y) = [J{y)m{y)] (2ciy + C2) = 2ciy + C2 . 
Then m{y) = 2c\y + C2 and = 65 = 0, so that 

ajii -I- 9712 -I- ajts -I- -I- ajis -I- mie = . 

10. 2"'' Elliptic Umbilic (5): Its corresponding 3-parameter induced map fc, for a given target point (si, S2) lying 
in the six-image region, is 

fc(a;, y) = {-2xy , -x^ + 5/ - 4ci/ - 3c2?/^ - 2c3y) = (si, S2) • 

Eliminating x, we obtain the polynomial 

^(y) EE -si - 4s22/' - 8C3/ - 12c2/ - IQciy'' + 20/ . 

The magnification is S[)T(a;, y) = 1/ {—Ax^ + 4,C3y + 12c2y^ + 24,cxy^ — AQy^). Substituting for via the equations 
ic{x,y) = (si,S2) gives 

9«(/= ^ ^ 



J{y) As2 + 12c3y + 24c2y2 + 40ci/ - 60/ 
One can check directly that v'(t/) = J{y){—2y), so that 

^'{y)m{y) = [J{y)m{y)] {-2y) = -2y . 
Then m{y) = —2y and 6„_i = 65 = 0, so that 

mi + dn2 + Ti3 + Dn4 + Ti5 + Dne = o . 
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11. 2"'' Hyperbolic Umbilic (5): Its corresponding 3-parameter induced map fc, for a given target point (si,S2) 
lying in the six-image region, is 

fc(a:, y) = {-2xy, -x^ - Sy" - 4ciy^ - 3c2?;^ - 2c3j/) = (si , S2) ■ 

Eliminating a;, we obtain the polynomial 

(^(y) = -s\ - 4s22/' - 8c3y=^ - I2C2/ - IGciy^ - 20/ . 

The magnification is 97l(a;, ?/) = l/(— 4a;^ + 4c3?/ + 12c2y^ + 24ci2/'^ + 40j/^). Substituting for via the equations 
fc(a;,y) = (si,S2) gives 

artfy) = ^— = ^- • 

J{y) 4s2 + 12c3y + 24c2y2 + 40ciy3 + 60?/4 

Once again, it is easy to check that tf'{y) = J{y){—2y), so that 

^'{y)m{y) = [J{y)m{y)] {-2y) = -2y . 
Then m{y) = —2y and &„_i = 65 = 0, so that 

+ aji2 + 9JI3 + + 9JI5 + arie = . 

This completes the proof. □ 

VI. CONCLUSION 

The paper presented a theorem about the magnification pre-images for caustic singularities up to codimension five. 
We proved that for generic families of general mappings between planes locally exhibiting such singularities, and for 
any point in the target lying in the region giving rise to the maximum number of real pre-images, the total signed 
magnification of the pre-images sums to zero. The signed magnifications are invariants as they are Gaussian curvatures 
at critical points. Our result extends earlier work that considered the case of singularities through to codimension 
three. The proof of the theorem is algebraic in nature and utilizes the Euler trace formula. In fact, we established a 
proposition that relates the coefhcients of the unique polynomial in the coset of certain rational funtions to Newton 
sums. It was then shown that the Euler trace formula follows readily as a corollary of our proposition. The findings of 
the paper are expected to be relevant to the study of dark matter substructures on galactic scales using gravitational 
lensing. In addition, since the results hold for generic general mappings, they are applicable to any system in which 
stable caustic singularities appear. 
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APPENDIX A: PROOF OF PROPOSITION [2] 

For convenience, we restate the result: 

Proposition. Consider any polynomial ip{x) = a^x" -I- • • • -I- aix -I- ao G C[a;] with distinct roots and any rational 
function h{x) G R. Let 

r{x) = bn-ix"^^ -I \-bix + bo 



be the respective unique polynomial representative of the coset ip'{x)h{x) in R/{ip{x)). Then the coefficients of r(x) 
are given in terms of the coefficients of and ifiix) through the following recursive relation: 

6„_i = c„_i6„„j^„_i H \- cibn-i,i + cobn-i,o i = l,...,n, (Al) 
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with 

bn-i,o = {n- {i- 1)) a„_(i_i) , i = l,...,n 



bn-i,k = bn-i,k-i + , i = 1, . . . ,n , k = 1, . . . ,n - 1 , 

an ^ 

where b-i^k-i = 0. 
Proof of Proposition. 

Wc begin with some preliminaries about quotient rings to make the proof more self-contained. Let C[x] be the ring 
of polynomials over C and let C(x) be the field of rational functions formed from quotients of polynomials in C[x]. 
The n zeros xi, . . . ,Xn of (p{x) = a„x" + ■ • • + aix + G C[x] are assumed to be distinct (generically, the roots of 
a polynomial are distinct). Let (f{x)) denote the ideal in C[x] generated by (p(x), and consider the quotient ring 
C[x]/{tp{x)), whose cosets we denote by g{x). This quotient ring has two important properties: 



• Property 1: If gi{x) = .92(2;), then by definition gi{x) — ,92(2;) = h(x)ip{x) for some h{x) G C[x], from which it 
follows that gi{xi) = g2{xi) for all n roots Xi of 1^3(2;). Thus members of the same coset must agree on the roots 
of ip{x), so that, in particular, J27=i9i-{xi) = J2"=i92{xi). 



Property 2: Each coset g{x) has a unique representative of degree at most n — 1, as follows: by the division 
algorithm in C[x], there exist polynomials q{x) and r{x) such that 

g{x) = q{x)^p{x) + r{x) , 



where deg r < deg (p = n. Passing to the quotient ring C[x]/{ip{x)), we see that g{x) = r{x). Suppose now that 
there exists another polynomial p{x) of degree less than n with g{x) = p{x). Then p{x) = r{x), so that 

p{x) — r{x) = h{x)ip{x) 

for some h{x) € C[a;]. If h{x) ^ 0, then deg hep > n, while the degree of the left-hand side is less than n. We 

must therefore have h{x) = and p{x) = r{x). We may thus represent every coset by its unique polynomial 
representative of degree less than n, which in turn implies that C[x]/{ip{x)) is a vector space of dimension n, 

with basis ^l,x,x^..., 

The next result will be used to show that Properties 1 and 2 also hold for a certain subset of rational functions in 
C{x) (see Claim 2 below). 

Claim 1. Let x\,...,Xn G C be distinct. Let ci, . . . ,c„ G C, not necessarily distinct. Then there exists a unique 
polynomial II{x) G C[x] with deg h < n such that H{xi) = Ci. 

Proof (Claim 1). Induction on n. For n — 1, define H{x) = c\. Now assume that the result is true for n — 1, 
and consider a set of n distinct complex numbers xi, . . . , a;„. By the induction hypothesis, there exists a polynomial 
h{x) G C[a;] with deg h < n — 1 such that h{xi) = c, for i = 1, . . . , n — 1. Now define 

H{x) = h{x) + ^ {Cn - h{Xn)) • 

[Xn - Xi){Xn - X2) ■ ■ ■ (Xn - Xn-l) 

It follows that H{x) G C[x] has degree less than n, and H{xi) = Ci for alH = 1, . . . , n. (As a simple example to show 
that H{x) need not be unique if the Xi, . . . ,Xn are not distinct, consider the numbers 2, 2, 3, 3 all being mapped to 0. 
Then the polynomials Hi{x) = {x - 2f{x - 3), H2{x) = (x - 2){x - 3f, and H^ix) = (x - 2)(x - 3) aU satisfy the 
assumptions of the lemma.) Suppose that there exist two polynomials Hi{x) and H2{x) with Hi{xi) = Ct = H2{xi). 
By the division algorithm in C[x], there are unique polynomials q{x) and r(x) such that 

iJi(x) - H2{x) = q{x) [(x - xi)(x - X2) • • • (x - x„)] + r(x) , 

where deg r < n. If q{x) ^ 0, then the degree of the polynomial on the right-hand side is at least n, whereas 
ili(x) — H2{x) has degree less than n. We must therefore have q{x) = 0. Moreover, if r{x) ^ 0, then Hi{xi) = H2{xi) 
gives that r{xi) = for all xi, . . . , x„. This implies, however, that r(x) has n distinct zeros and so must have degree 
n, a contradiction. Thus ifi(x) = H2{x). □ (Claim 1) 
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Let R C C(a;) denote the subring of rational functions that are defined at the roots Xi of 1^5(2;), 

1 pixj 

i? = < : p{x), q{x) G C[a;] and q{xi) ^ for all roots Xi 

I ^\-^) 

and consider the quotient ring R/{(p{x)). The next claim states that the ring R/{(p{x)) satisfies Properties 1 and 2. 

Claim 2. Members of the same coset in R/{(p{x)) agree on the roots Xi of (p{x), that is, if 51(2;) and g2{x) belong to 
the same coset, then gi{xi) — g2{xi), and so X]r=i5i(*^«) ~ X^ILi ffsl^^i)- In addition, any rational function h{x) G R 
will have in its coset h{x) G R/ {ip{x)) a unique polynomial representative r{x) of degree less than n. 



Proof ( Claim 2). Notice that, if hi{x) — h2{x) G R/ {ip{x)), then by definition there exists a rational function h{x) G R 
such that 

hi{x) ~ h2{x) = h{x)tp{x) , 

so that hi{xi) — h2{xi) for all the zeros Xi of >p{x). In other words, R/{ip{x)) also satisfies Property 1. It turns out 
that when the zeros xi, . . . ,Xn of f{x) are distinct, as we are assuming they are, then R/ {ip{x)) also satisfies Property 
2 (in fact R/{ip{x)) and C[x]/ {(p{x)) will be isomorphic as rings). For given a coset h{x) G R/{(p{x)), Claim 1 shows 
that there is a unique polynomial g{x) G C[a;] of degree less than n whose values at the n roots Xi are h{xi). Then 
the rational function g{x) — h{x) G R vanishes at every Xi, and a simple application of the division algorithm applied 
to the numerator of g{x) — h{x) shows that g{x) — h(x) G R/{ip{x)). Thus any rational function h{x) G R will have 
in its coset h{x) G R/{if{x)) a unique polynomial representative r{x) of degree less than n. □ (Claim 2) 

We now begin the proof of the Proposition by establishing the following Lemma: 

Lemma. Let tf{x) = a„a;" + ■ • ■ + aix + oq and consider the quotient ring R/ {Lp{x)). For any 1 < fc < n — 1, let 

rk{x) = bn-i,k a;""^ H h foi,fe x + bo,k 



be the unique polynomial representative in the coset Lp'{x)x^ . Then the following recursive relation holds: 
K-ifi = {n- (i- 1)) a„_(,_i) , i = 1, . . . , n , 

^n-i,fc = &n-i,fc-i + 1—1 , i = l,...,n, fc=l,...,n-l, 

an 

where b-i^k-i = 0. 

Proof of Lemma. The existence and uniqueness of the polynomial 

rkix) = fe„_i,fcx""^ H h 6„_i,fcx""* H h biMX + bo,k 1 , 

where 



(A3) 



ip'{x)x'' ^rk{x) ^bn-ukx"-^ + --- + bn-t,kx"-' + --- + bi.kX + bo.k'^ , (A4) 



were established in Claim 2. Also, note that since ip{x) = G R/{ip{x)), we have 



^x-^T. (A5) 



Case fc = 0: By (gU, we get 



Lp'{x)x° = rQ{x) = bn-i.ox'^ h &„-i,o a:" ' H h 61,0 a^ + 60,0 1 

However, 



ip'{x)x^ — f'{x) = na„ x" ^ + • • • + (n — (i — l))a„_(j_i) x" * + ••• + 2a2 a; + ai 1 . 
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Consequently, 



bn-i,o = {n- [i- 1)) a„_(i_i) 
Case /c = l,...,n— 1: Equations (|A4|) and (|A5|) yield 



ip'{x)x'' = fe«-i,fca;" h 6„-i,fc h 6i,fe .x + 5o^fc 1 



fe-i 



= a; 6„_i,fe_i.x" 1 + 6„-2,fc-i a;" h 6i,fe_i x + 6o^a;-i 1 

= bn~i.k-i + bn-2,k-i a;""^ H h foi,fe-i a::^ + bo^k-i x 

— bn-\,k-\ 



an 



ai _ flQ _ 

X 1 



bn-2.k-lX'^ ^ 



E 



■ bn^l,k~l + &n-(i+l),fc-l 



The coefficients of (jA4l) are then related to the coefficients of ai of 'p{x) as follows: 



(A6) 



x"^ + 6o,fc-l X 



^n—i.k — 



■ 6„_i,fc_i + 6„_(i+i)^fe_i , «=l,...,n, A: = l,...,n-1, 



where the coeffiencients b^-ifl are given by (jA6p . Note that 6„ ^ = since the unique polynomial goes up to degree 
n — 1. □ (Lemma) 

We now complete the proof of the Proposition. If hi^^{x) and h2.*{x) are the unique polynomial representatives 
of the cosets hi{x) and h2{x), respectively, then by uniqueness, the sum /ii^*(a;) + h2,*{x) is the unique polynomial 
representative of the coset hi{x) + h2{x). With that said, we note that, since h{x) = h^,{x), it follows that r{x) = 
ip'{x)h{x) = (p'{x)h^{x). We thus have 



r{x) — tp'{x)h^,{x) 



c„_i</5'(x)a:" 1 + • • • + ciip'{x)x + co(p'{x) 
Cn-ir„-i{x) H h ciri{x) + coro(x) 

n n 

Cn-i ^ 6„-i,„-ia;"^* H ^ ci ^ 6„- 



(c„_i6„_i,„_i H h ci6„_,;4 + co6„-i,o) a:" 

n 

^^6„_ia;"~* . □ (Proposition) 
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